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Abstract. Using the University of Hawaii's 8K mosaic camera (UH8K), we have measured the angular correlation 
function lj(8) for 100,000 galaxies distributed over four widely separated fields totalling ~ 1 deg 2 and reaching 
a limiting magnitude of Iab{3o, 3") ~ 25.5. This unique combination of areal coverage and depth allows us to 
investigate the dependence of u)(6) at 1', j4 w (1'), on sample median magnitude in the range 19.5 < lAB-med < 24. 
Furthermore, our rigorous control of systematic photometric and astrometric errors means that fainter than 
lAB-med ~ 22 we measure u)(d) on scales of several arc-minutes to an accuracy of 30%. Our results show that 
A u (l') decreases monotonically to Iab ~ 25. At bright magnitudes, lu{9) is consistent with a power-law of slope 
8 = —0.8 for 0.2' < < 3.0' but at fainter magnitudes we detect a slope flattening with 5 ~ —0.6. At the 3<r 
level, our observations are still consistent with 5 = —0.8. We also find a clear dependence of Au,(l') on observed 
(V — I)ab colour. In the magnitude ranges 18.5 < Iab < 24.0 and 18.5 < Iab < 23.0 we find galaxies with 
2.6 < (V — I)ab < 2.9 (the reddest bin we consider) have A u (l')'s which are ~ lOx higher than the full field 
population. On the basis of their similar colours and clustering properties, we tentatively identify these objects 
as a superset of the "extremely red objects" found through optical-infrared selection. We demonstrate that our 
model predictions for the redshift distribution for the faint galaxy population are in good agreement with current 
spectroscopic observations. Using these predictions, we find that for low-fi cosmologies and assuming a local 
galaxy correlation length ro — 4.3/t -1 Mpc, in the range 19.5 < lAB-med < 22, the growth of galaxy clustering 
(parameterised by e), is e ~ 0. However, at 22 < lAB-med < 24.0, our observations are consistent with e > 1. 
Models with e ~ cannot simultaneously match both bright and faint measurements of Au,(l'). We show how 
this result is a natural consequence of the "bias-free" nature of the "e" formalism and is consistent with the field 
galaxy population in the range 22.0 < Iab < 24.0 being dominated by galaxies of low intrinsic luminosity. 
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1. Introduction 

Understanding the formation of structure in the Universe 
is one of the most pressing questions in modern cosmology. 



The Sloan and 2dF surveys currently in progress (Colless 



1998; Gunn 1995) will provide an accurate picture of large- 



scale structures in the local Universe but presently our 
knowledge of galaxy clustering at z > 0.5 is poorly con- 
strained. This is entirely a consequence of the technical 
limitations in imaging and spectroscopic equipment, which 
(until recently) have had fields of view ~ 50 arcmin 2 ; in 
all cosmologies this translates to < 1 h~ 1 Mpc at z > 0.5. 
Covering a substantial enough area to provide meaningful 
statistics on larger (10 — 20/i~ 1 Mpc) scales at higher red- 



H. J. McCracken et al.: The Canada- France deep fields survey-I 



shifts (z ~ 1) has been prohibitively expensive in telescope 
time. Consequently, many galaxy clustering measurements 
made at these redshifts have been dominated by the ef- 
fects of sample variance, and also have only been able to 
investigate the highly non-linear regime where the predic- 
tions of theoretical models for the clustering of galaxies 
are strongly dependent on the biasing schemes employed. 
Other studies, such as investigating the variation of clus- 
tering amplitude by galaxy type or intrinsic luminosity, or 
the accurate measurement of higher-order statistics such 
as S3 have been even more challenging. 

However, with the advent of wide-fiel d multi-object 



spec trographs like VIRMO S and DEIMOS ( |Le Fevre et al 



2000; Davis fc Faber 1998), in addition to wide-field mo- 



saic cameras on 4m-telescopes, this situation is changing. 
In this paper we detail a new survey, the Canada-France 
Deep Fields (CFDF) project which has been carried out 
using the University of Hawaii's wide field (28' x 28') 8K 
mosaic camera, UH8K. This survey has targeted three of 
t he original field s of the Canada-France redshift survey 
( Lilly et al. 1995 ). In total the CFDF consists of four in- 
dependent deep fields, each of area 0.25 deg 2 . All of these 
have VI colours, three BVI and two and a half UBVI. 
The survey reaches a limiting magnitude (3ct, 3" aperture) 
of Iab ~ 25.5 and at least one magnitude fainter in UBV 
(Table W. The ~ 10 5 galaxies in the survey, coupled with 
1000 spectroscopic redshifts present throughout our fields, 
allows us to investigate with unprecedented accuracy the 
evolution of galaxy clustering to z ~ 1 (the survey's me- 
dian redshift at its completeness limit of Iab ~ 25.5). 
Moreover, our four widely separated fields also ensure that 
we can estimate the effect of cosmic variance on our re- 
sults. 

To date, there have been many studies of U){6) carried 
out using deep imaging surveys conduc ted using charge- 



cou pled-device (CCD)-based dete c tors ( McCracken et al 
2000|; |Woods fe Fahlman 19971; |Hudon fe Lilly 1996| ; 



Brainerd et al. 1994 ; Roche et al. 1993 ) . These works have 
generally focussed on one or two fields, usually covering 
~ 50 arcmin 2 each and typically reaching limiting magni- 
tudes of Iab ~ 25. Several authors have also attempted 



tions. One can then conclude which value of e is most 
appropriate for any given set of observations. Based on 
comparisons between uj{9) measurements in deeper CCD 
surveys and photographic measurements at brighter mag- 
nitudes, many authors concluded that, for z < 1 at least, 
growth of galaxy clustering was consistent with < e < 2 
(Brainerd et al. 1994). More recently, direct measurements 
of ro(z) have been attempted at z < 1 using spectro- 
scopic samples (Carlberg et al 



2000t Small et al. 1999 



Le Fevre et al. 1996; Cole et al. 1994). These works have 



1 



demonstrated the importance of sample selection in mea- 
suring galaxy clustering evolution; e has been shown to be 
sensitive to the range of intrinsic galaxy luminosities and 
spectral types selected. Attempts have also been made us- 
ing photometric redshifts computed using either ground- 
based or space -based imaging data to measure the growth 
of clustering (|Tcplitz et al. 200l|; |Brunner et al. 200"0; 



Arnouts et al. 1999a ; Connolly et al. 1998 ). However, 
the finding that clustering amplitudes for Lyman-break 
galaxies was similar to some classes of galaxies found lo- 
cally ( Adelberger et al. 1998| ; Giavalisco et al. 1995 ) has 
provided the clearest evidence to date that this simple 
formalism could not fully account for the observations of 
clustering at z ~ 3. 

In this paper, the first in a series, we will introduce 
the CFDF survey, explain in detail our data reduction 
strategy and demonstrate its robustness. As a first appli- 
cation of this dataset, we will present a measurement of 
the projected galaxy correlation function uj{9). The angu- 
lar size and depth of the CFDF allows us to make a reli- 
able determination of u>{0) over a large magnitude range 
(18 < lAB-med < 24). Moreover the four separate fields 
allows us to make an estimate of the field-to-field variance 
in the galaxy clustering signal. Finally we will discuss how 
appropriate the "epsilon" formalism is to describe the evo- 
lution of galaxy clustering measured in our data. 

In a future paper (Foucaud et al., in preparation) we 
will describe our measurements of the clustering length tq 
at z ~ 3 from a sample of ~ 2000 Lyman-break galaxies 
derived from the CFDF dataset. By adding R— and Z— 



band data from the new CFH12K camera (Starr et al 



many sepa rate pointings ( Roche fe Eales 1999 ; Postman 
et al. l |998] ), although these surveys reach much shallower 
limiting magnitudes {Iab ~ 22). In contrast, the CFDF 
survey, by virtue of its depth and angular coverage, is able 
to provide an accurate measurement of ui{6) in the range 
18.5 < Iab < 25.0. 

Normally the results from these surveys have been in- 
terpreted in terms of the "e— " formalism first introduced 
to explain clustering amplitudes observed at bright mag- 



to cov er larger areas (> 1 deg^) by mosaicing together 200C ) we expect to sufficiently increase the accuracy of 



nitudcs on photographic plates (Groth & Peebles 1977; 



Phillipps et al. 1978). With this approach, an assumed 
redshift distribution dN/dz (or one measured from an in- 
dependent spectroscopic survey) and cosmology is coupled 
with a model for the evolution of £(r, z) (parametrised 
by e). In this way it is possible to predict the amplitude 
of u)(0) at any magnitude limit, based on these assump- 



thc photometric redshifts in the range < z < 1 to allow 
a direct measurement of tq(z) in this interval; however, in 
this paper we will concern ourselves only with measure- 
ment of u>(6) and its dependence on apparent magnitude 
and colour. 



2. Observations and reductions 

2.1. Observations 

B, V , I observations were taken on the Canada-Francc- 



Hawaii teles cope with the UH8K mosaic camera ( Mctzgci 
|ct al. 1995 ) over a series of runs from December 96 



- June 97; details are given in Table [y. Typically, for 
the VI— band exposures we used exposures of 1800s; 
for the B— images we adopted exposure times of 2400s. 
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Individual exposures which had a FWHM > 1.2" were dis- 
carded. Observing conditions were generally quite stable: 
for example, for the 03hr field observations In V and / 
bands, the median seeing is ~ 1.1" and 1.2" respectively. 
Additionally, as our point-spread function (PSF) is almost 
always oversampled, it is not necessary to carry out PSF 
homogenisation before image stacking. At each pointing 
there is ~ 10 exposures which allows us to carry out ad- 
equate cosmic ray removal and to fill the gaps between 
each CCD in the mosaic. 

The UH8K camera consists of eight frontside- 
illuminated Loral-3 2048 x 4096 CCDs, arranged in two 



high and low darks are identified by their statistics (mean, 
median). Two full reductions are done. Each kind of dark 
is subtracted and then a dark-independent flat (dome or 
twilight) is applied. The flatter final image (with small- 
est amplitude of residual flatness variations) indicates the 
kind of dark that was actually present in the data. 

We generate "superflats" from the science images 
themselves as dome flats or twilight flats by themselves 
produce residual sky variations > 1%. These superflats 
are constructed by an iterative process, which begins by 
the division of our images by a twilight flat. On these 



banks )f four devices each. Each bank is read out sequen- 
tially. The upper-right CCD (number 8) has very poor 
charge-transfer properties and data from this detector was 
discarded. The pixel scale is 0.205"pixel _1 . Additionally, 
all of the CCDs have separate amplifier and controller elec- 
tronics. This arrangement, in addition to the necessity of 
removing the CFHT-prime focus optical distortion before 
stacking our images, resulted in a lengthy data reduction 
procedure which is outlined in the following sections. 

Because of the poor blue response of the Loral-3 de- 
vices, U— band observations for the CFDF survey were 
taken with the Kitt Peak 4m Mayall telescope and the 
Cerro Tololo Inter- American Observatory's (CTIO) 4.0m 
Blanco Telescope during a series of runs in 1997. Because 
of the smaller field of view of these cameras, four separate 
pointings were needed to cover each UH8K field. While 
the seeing on the BVI frames is ~ 0.7" — 1.1" some of 
the U— stacks are significantly worse (1.2" — 1.4") which 
had to be accounted for during catalogue preparation (two 
catalogues were prepared: one for those science objectives 
which required U— data and one for those which did not; 
this is explained in more detail in Foucaud et al.). 

2.2. Preprocessing 

Pre-processing followed the normal steps of overscan cor- 
rection, bias subtraction and dark subtraction. We fit a 
fourth-order Legendre polynomial to the overscan region 
to allow us to remove structure in the overscan pattern. 
The high dark current (~ O.le^s -1 ) of the UH8K makes it 
essential to take dark frames. For I— and V— bandpasses, 
where the sky background is high, one unique dark frame 
(composed of the average of 5-10 individual exposures) 
can be used; however, with the B— band special steps 
have to be taken due to the '"dark-current jumping" ef- 
fect which manifests itself as a either "high" or " low" dark 
current level, which can affect either right or left banks 
of CCDs independently. Because of the low quantum ef- 
ficiency of the UH8K CCDs in B the dark current is a 
significant fraction of the sky level and consequently ac- 
curate dark subtraction must be performed to produce 
acceptable results. We achieve this by generating two sets 
of darks: "high" darks and "low" darks which we apply 
by a trial-and-error method to each B— exposure to deter- 
mine which dark is appropriate for a given dataset. The 



twilight-flattened images we run the sextractor (Bertin 
fc Arnouts 1996]) package to produce mask files which 



identify the bright objects on each frame. For large sat- 
urated stellar objects we grow these masks well into the 
wings of the point-spread-function by placing down circles 
on these objects. In addition, we mask out non-circular 
transients (typically scattered light and saturated columns 
near bright stars) on some images. Using these mask files 
we combine each image to produce a superflat each pixel 
of which contains only contributions from the sky and not 
object pixels. After the division of the twilight-flattened 
images by the superflats, the residual variation in the sky 
level is < 1%. 

Because the gain and response of each CCD in the 
mosaic is not the same, we scale our flat fields for each 
filter so that the sky background in each chip after division 
by the flat field is the same ( norm ally these exposures 
are scaled to chip 7). In Sect. 2.5 we will quantify how 
successful this procedure is in restoring a uniform zero- 
point over the entire field of view of the image. 

2.3. Astrometric image mapping 

For each set of observations in each filter of our field, we 
have typically ~ 10 pointings (we use the term "pointing" 
to refer to the eight separate images which comprise each 
read-out of the UH8K camera). Each of these pointings 
are offset by ~ 5" from the previous one in a random 
manner; these offsets allows us to remove transient events 
and cosmetic defects from the final stacks, and also to 
ensure that the gaps between the CCDs (which are ~ 3") 
are fully sampled . Additionally, on some of our fields, 
pointings were taken over several runs with the camera 
bonette in different orientations. Given the non-negligible 
optical distortion at the CFHT prime focus (amounting 
to a displacement of several pixels at the edge of the field 
relative to an uniform pixel scale) this means it is essential 
that these distortions are removed before the pointings can 
be coadded to produce a final stack. A further requirement 
is that each of the stacks for each of the filters can also 
be accurately co-aligned, for the purposes of measuring 
aperture colours reliably. 

In the mapping process, the images from each CCD are 
projected onto an undistorted, uniform pixel plane. The 
tangent point in this plane is defined as the optical centre 
of the camera, and is the same for each of the eight CCDs. 
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Field R.A. (2000) Dec. (2000) Band Exposure time Seeing 3<r measurement Area Date 

(hours) (arcsec) (AB mags) (deg 2 ) 



0300+00 


03:02:40 


+00:10:21 


U 


10.8 


1.0 


26.98 


0.25 


06/97 








B 


5.5 


1.1 


26.38 


0.25 


09/97 








V 


4.2 


1.3 


26.40 


0.25 


12/96 








I 


5.5 


1.0 


25.62 


0.25 


12/96 


2215+00 


22:17:48 


+00:17:13 


u 


12.0 


1.4 


27.56 


0.12 


06/97 








B 


5.5 


0.8 


25.90 


0.25 


09/97 








V 


2.7 
1.5 


1.0 


26.31 


0.25 


06/97 
09/97 








I 


1.7 
2.7 


0.7 


25.50 


0.25 


06/97 
09/97 


1415+52 


14:17:54 


+52:30:31 


u 


10.0 


1.4 


27.71 


0.25 


03/97 








B 


5.3 


0.8 


26.23 


0.25 


06/97 








V 


2.3 


1.0 


25.98 


0.25 


06/97 








I 


2.6 


0.7 


25.16 


0.25 


06/97 


1130+00 


11:30:02 


-00:00:05 


V 


3.3 


0.8 


26.42 


0.25 


05/97 








I 


4.4 
3.0 


1.0 


25.80 


0.25 


12/96 
05/97 



Table 1. Details of the images used in the CFDF fields. For each field we list the total integration time, in addition to 
the 3cr detection limit inside an aperture of 3". For the 22hr and llhr fields, the fields are composed of two separate 
stacks with bonette rotations, as detailed in the text. For the 22hr field the U— data covers only half the field. 



Overall, our goal is to produce an root-mean-square reg- 
istration error between pointings in each dither sequence 
and between stacks constructed in different filters which 
does not exceed one pixel (0.205") over the entire field of 
view. 

Our astrometric mapping process is essentially a two- 
step process. In outline, this involves first using the United 
States Naval Observatory (USNO) catalogue (Monet] 
1998y~Jo derive an absolute transformation between (x,y) 
(pixels) to (a, S) (celestial co-ordinates). Following this, 
a second solution is computed using sources within each 
field. This method ensures that our pointings are tied to 
an external reference frame, but also provides sufficient 
accuracy to ensure that pointings can be registered with 
the precision we require; the surface density and positional 
accuracy of the USNO stars is too low to ensure this. To 
fully characterise the distortions which are present in the 
camera optics we adopt a higher order-solution, which con- 
sists of a combination of a standard tangent plane projec- 
tion and higher-order polynomial terms. To prevent so- 
lution instabilities at the detector edges we use a third- 
order polynomial solution. To compute the astrometric 
solutions and carry out the image mapping we use the 
mscred package provided within the IRAF[J data reduc- 
tion environment. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 



For each field we begin our procedure with the I— 
band, as these exposures normally have the highest num- 
bers of objects. Using the external catalogue we com- 
pute an astrometric solution containing a common tan- 
gent point for each of the eight CCDs. Typically, we find 
50 - 100 sources per CCD, with a fit RMS of - 0.3". Next, 
using the task mscimage we project each of the eight im- 
ages onto the undistorted tangent plane, using a third- 
order polynomial interpolation. Following this, we extract 
the positions in celestial co-ordinates of a large number 
(~ 1000) of sources distributed over all eight CCDs. This 
list forms our co-ordinate reference, and we use this list 
in conjunction with the task mscimatch to correct for the 
adjustments in the WCS (world co-ordinate system) due 
to slight rotation and scale change effects for each succes- 
sive pointing in the dither. Before stacking we also remove 
residual gradients by fitting a linear surface and scale the 
images to photometric observations if necessary. Because 
each image now has a uniform pixel scale we need only 
apply linear offsets before constructing the final stack. 
Setting the gaps between each CCD to large negative val- 
ues which are rejected in the stacking process allows the 
production of a final, contiguous image. To stack our im- 
ages we use a clipped median, which although not optimal 
in signal-to-noise terms, provides the best rejection of out- 
lying pixels for small numbers of pointings. 

From this final, combined stack, we extract a second 
catalogue of (~ 1000) sources (with a, S computed from 
the astrometric solution) which we use as an input 'as- 
trometric catalogue' for the dither sequences observed in 



H. J. McCracken et al.: The Canada-France deep fields survey-I 

Field 6 



0.5 






0.5 



-1 



T 



x ■. 



■ -- ■ ■ i 

• ■" ■<■ -ifs ■-*■:■* ".". 

■ ■ ■ ^ * *? .p J 

■■■ -fA." \ k '. . 
■■ ^\. ■ ■ . 



j_ 



-0.5 




<5x 



0.5 



Fig. 1. The difference, in pixels, of the positions of non- 
saturated stellar sources with 18 < / < 22 in the 03hr I— 
band stack compared with the 03hr B— stack. The rms in 
both co-ordinate directions is ~ 0.3 pixels, or ~ 0.06". 



other filters (as opposed to the USNO catalogues which 
we use for the first step). Typically, the RMS of the fit 
these cases is < 0.1". We then proceed as before, map- 
ping each of CCD images from each pointing in the dither 
set onto the undistorted tangent plane and constructing a 
final stack. 

For the final mapping between the stacks taken in dif- 
ferent filters, we find a residual of ~ 0.06", or ~ 0.3 pixels 
over the whole field of view, which is with our aim of a 
root-mean-square of one pixel or less; this is illustrated in 
Fig. Q. Our large grid of reference stars extracted from the 
I— band image ensures that the derived WCS for the other 
filters is very well matched to the I— band exposure. We 
also find that this method allows us to successfully reg- 
ister and combine observations distributed over separate 
runs containing bonette rotations. 

For the U— band exposures, each of the four corners 
were stacked separately and scaled to have the same pho- 
tometric zero-point. Then, using the I— band reference 
list described previously, a mapping was computed be- 
tween each of the corners and the undistorted I— stack. 
During this process the image was also resampled (using 
the same third-order polynomial kernel employed above) 
to have the same pixel scale as the UH8K data. These 
four U— images were then stacked to produce the final 
U— mosaic. Overall, we find that the rms of the mapping 
between I— and U— is not as good as between the UH8K 
data, with an rms ~ 1 pixel in the region of CCD 8 (the 
CCD suffering cosmetic defects) but still within our stated 
goal. 



I E{B - V) a E(B - V) b 

Schlegel et al. BH 



0300+00 -48 179 

2215+00 -44 63 

1415+52 +60 97 

1130+00 +57 264 

Table 2. The (I, b) for each CFDF field, together with the 



0.071 


0.040 


0.061 


0.040 


0.011 


0.000 


0.026 


0.013 



galactic dust extinction corrections from Schlegel et al 



(1998)( a ) compared to the values of Burstcin & Heiles 
(1982) ( 6 , BH). 



Each of the final images have a scale of 0.204" /pixel 
and cover ~ 8000 x 8000 pixels (the scale of the final 
stack is determined from the linear part of the astrometric 
solution of the image which is closest to the tangent point 
of the camera). In all the analyses that follow we exclude 
the region covered by CCD 8 as this chip has very bad 
charge transfer properties and is highly photometrically 
non-linear. However, for the llhr-I and 22hr-I stacks we 
are able to use the full area of UH8K because these final 
stacks consist of two separate stacks with bonette rotation, 
allowing us to cover the region lost by the bad CCD. 

2.4. Calibration 

In this Section we will describe how we derive the relation- 
ship between magnitudes measured in our detector/filter 
combination (which we denote by u c fdf, hfdf, v c fdf,i c fdf) 
and the standard Johnson UBVI system. Our zero-points 
are computed from observations of the standard star fields 
of Landolt (1992Q . Of the four observing runs with UH8K 
which are discussed here, only the data from May were not 
photometric and for the runs of June and October suffi- 
ciently large numbers of observations of standards were 
taken (n ~ 30) it was possible to determine the colour 
equation. 

We apply the same data reduction procedure to the 
standard star fields as we do to the science frames. This in- 
volves bias and dark subtraction, followed by flat-fielding 
which is necessary to account for the sensitivity and gain 
variations from CCD to CCD and the application of the 
astrometric solution derived previously to produce a sin- 
gle, undistorted image. This procedure assures that a uni- 
form pixel scale is restored when computing the photomet- 
ric zero-point, and has added advantage of that we may 
use a catalogue of standards in a, S to derive the zero-point 
in a semi-automated fashion. All standards used were visu- 
ally inspected and faint or saturated objects were rejected. 
Our zero-points are corrected for galactic extinction using 
the E{B - V) values provided by jBchlegel et al. (1998| ). 

In Fig. g we show sample plots of standard star obser- 
vations taken during the October observations. For the 
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Fig. 2. Standard star calibration plots for the I and V 
filters for the October 97 data (upper and lower panels 
respectively). For each panel, we plot the difference be- 
tween the standard star magnitude and the instrumental 
magnitude (corrected to one airmass) as a function of the 
true I, V magnitudes. The solid line shows the adopted 
zero-point; this is 24.24 ± 0.04 (V) and 24.26 ± 0.06 (J), 
per second and at one airmass. Atmospheric extinction co- 
efficients are taken from the CFHT handbook, although 
most of our standards are at or near the zenith. 



I and V band for all three runs we derive zero-point 
rms of ~ 0.05 magnitudes, with no evidence for position- 
dependent residuals (as would happen if an error had oc- 
curred in the flat- fielding process). 

Our observations do not indicate that the presence of 
a colour term for either the V or I filters and in what 
follows we assume V = v c fdf and I—i c fdf- However, we 
find that b c fdf is different from the standard Johnson B— 
and for this reason we derive colour terms. 



°cfdf 



B - 0.07(5 - V) - 23.2. 



(1) 



We convert our magnitudes to the AB system using 
Uab = u c fdf + 0.73, B A b = b cfd f -0.1, V AB = v cfd f and 
Iab — icfdf + 0.43 (computed based on our filter response 
functions). 

2.5. Systematic photometric errors 

To measure accurately the galaxy clustering signal it is es- 
sential that the photometric zero-point is uniform across 
the stacked images. Zero-point variation across the mosaic 
will introduce excess power on large scales and contribute 
to a flattening of the u>{6) on large scales. For single-CCD 
images, improperly flattened data can produce this effect; 



Fig. 3. B — b c fdf vs (B — V), based on observations of 
51 standard stars during the October 1997 observations. 
The slope and offset of the fitted line is 0.07 and 23.26, 
respectively. The slope is not well determined due to the 
small range in (B-V) colours spanned by our observations 
and we expect errors of ~ 0.02. 

in our case we have the additional complication that we 
must correctly account for the different responses and am- 
plifier gains for the eight CCDs in our mosaic. As outlined 
above, this is accomplished by scaling each CCD image 
before co-addition to have the same sky back ground. Our 
standard star reductions detailed in Sect. |]j have already 
indicated that this procedure produces zero-point varia- 
tions on order 0.05 magnitudes r.m.s. . However, further 
observations allow a more rigorous test of this effect. On 
two separate occasions we have observed the same field 
(llhr, 22hr) after the camera bonnette had been rotated 
180 deg. These data provide an excellent opportunity to 
verify that there are no residual systematic magnitude ze- 
ropoint variations in our final stacks after co-addition and 
stacking have been carried out. 

To carry out these tests we prepare two separate 
stacked mosaics. For the field at llhrs, we have 4.4hrs of 
integration in I from December 1996 and 3hrs total inte- 
gration from May 1997. By using sources extracted from 
the December run to compute our astrometric solution 
following methods outlined above, we can produce final 
stacked mosaics which are aligned with a standard devia- 
tion of < 1 pixel over the entire field of view. By carrying 
out the detection process on the sum of these two images, 
and photometry on the two separate mosaics, we are able 
to measure the difference in magnitude between sources 
located in the same part of the sky but falling on differ- 
ent elements of the detector-telescope system. Note that 
because we place our photometric apertures on the same 
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Fig. 4. The magnitude difference as a function of x and y 
position between non-saturated stellar sources with 18.5 < 
Iab < 22.5 in the December 1996 llhr stack and the May 
1997 data covering the same field. We find a(S m ) ~ 0.04 
magnitudes. 



positions on each of the two stacks, this test also allows 
us to investigate magnitude errors introduced by mapping 
inaccuracies between the two stacks, which are expected 
to be present for the measurement of aperture colours. 

The results of this test are illustrated in Fig. |] where 
we plot the difference in magnitude for non-saturated stel- 
lar sources with 18.5 < Iab < 22.5 between the two llhr 
stacks as a function of position in both x and y directions. 
We find that the systematic magnitude errors, measured 
as the dispersion of these residuals is ~ 0.04 magnitudes, 
which corresponds to the limit of our CCD-to-CCD cali- 
brations, as explained in Sect. lL4 . 

Towards the 50% completeness limit of our catalogues, 
Iab ~ 25.5, differential incompleteness becomes a signif- 
icant bias in the measurement of u)(6). This effect arises 
from the differing read-out electronics and detector gains 
used in each of the eight individual CCDs in UH8K. 



Neuschaefer & Windhorst (1995), using the Palomar four 



shooter camera, discuss this effect in more detail. However, 
we emphasise that all our scientific analysis is only carried 
out where our completeness is > 80%, as determined from 
the simulations and source counts detailed in the following 
sections. Furthermore have verified that this effect is only 
significant at the faintest magnitudes by adding 40,000 
objects with the same clustering amplitude as galaxies at 
Iab ~ 25 to o ne of our images. This test is described in 
detail in Sect. 5.2. 



2.6. Random photometric errors 

We may also use these repeated observations of the same 
field to investigate what random photometric errors are 
present in our data. At fainter magnitudes these errors 
dominate. We have used three separate methods to esti- 
mate the magnitude errors computed in our data; firstly, 
we may use the errors computed directly by sextractor; 
secondly, we can use the errors computed from the sim- 
ulations detailed in Sect. 2.7 in which stellar objects are 
added to our fields and recovered; and lastly, we may use 
the our two independent stacks of the same field to esti- 
mate our errors. 

Fig. shows magnitude errors for these three different 
estimators: sextractor (circles), the simulation (stars) 
and from the direct measurement (squares). We have car- 
ried out these tests on both the llhr stacks and the 
22hr stacks. In many magnitude ranges, the sextractor 
errors are lower than the other two measurements. We 
believe the origin of this discrepancy is due to the im- 
age resampling and interpolation process which produces 
images with correlated background noise. By contrast, 
the sextractor magnitude errors are computed assum- 
ing white background noise. 

2.7. Incompleteness simulations and limiting 
magnitudes 

In Table [y we list the 3c values for detection in a 3" 
aperture. These should be regarded as lower limits on the 
detectability of the galaxies in our catalogues. To better 
characterise the photometric properties of our images we 
have carried out an extensive set of simulations. These 
simulations involved adding artificial stars and galaxies to 
our single-band images and measuring the fraction recov- 
ered as a function of magnitude. In Fig. O we show the 
results of one set of such simulations for the 03hr field. 

We note that this result should be regarded as lower 
limit to the completeness in our data as our actual cata- 
logues are constr ucte d using the chisquared technique de- 
scribed in Sect. 34 and can be expected to be slightly 
deeper (but note also that in constructing the chisquared 
catalogues all images must first be convolved to the worst 
seeing). From a rough comparison of the /—band galaxy 
counts presented in Fig. g, we see that the we see that the 
simulations provide a good estimate of the magnitude at 
which the the observed counts begin to fall off. 

2.8. Comparison with CFRS photometry 

Three of our fields (22hr, 03hr, 14hr) cover the original sur- 
vey fields of the C anada-France Redshift Survey (CFRS; 
Lilly ct al. (1995[ )). For the 22hr and 14hr field we have 



BVI photometry from the CFRS; for the 03hr field, data 
exists in the VI bandpasses. 

For all these fields we have carried out a detailed com- 
parison of our photometry with CFRS photometry. In 
Fig. |7J we compare V and / photometry from our stacked 
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Fig. 5. Root-mean-square (RMS) magnitude errors as a 
function of Iab magnitude. The upper and lower panels 
shows measurements made on the llhr and 22hr fields 
respectively. Magnitude errors were computed using a va- 
riety of techniques: from sextractor (open circles); from a 
simulation in which stars were added to random patches 
of the field (stars) and the dispersion in the recovered 
total magnitudes is calculated; and finally as measured 
between two independent stacks covering the same field 
(open squares). The direct magnitude error measurements 
have been multiplied by v2 to account for the shorter ex- 
posure time for the individual stacks. At all magnitudes 
the sextractor errors appear to underestimate the both 
the direct and simulation error by at least a factor of two. 



images with the CFRS for V and / filters in the 03hr 
and 14hr fields. For the 14hr fields, the agreement with 
the CFRS photometry is < 0.1 magnitudes or better. In 
the 03hr field, however, we find that our magnitudes are 
~ 0.2 and 0.1 magnitudes brighter than CFRS magnitudes 
in the V and / filters respectively. We suspect the origin 
of this discrepancy is that the CFRS fields were selected 
to have low galactic extinction as measured in the maps of 
Burstein fc Heiles (1982| ) (BH). In Table § we show that 
the difference between the BH extinction and the more 
recent E(B — V) values given in [Schlcgcl ct al. (199S| ) is 
non- negligible (amounting to ~ 0.15 in Iab magnitudes). 
In all our fields we ap ply extinction correct ions based on 
E(B-V) values from |5chlcgcl ct al. (1998J) . 



3. Catalogue preparation 

3.1. Construction of merged catalogues 

To prepare merged catalogues with UBVI colours for each 
object we derive a "detection" image to locate objects and 



Fig. 6. Derived completeness in four bands for the 03hr 
field as a function of AB magnitudes. These simulations 
were derived by adding 1,000 artificial stars in each half- 
magnitude interval to the original images and measuring 
how many objects were recovered as a function of magni- 
tude. The detection threshold used for this test (number of 
sigmas above the noise background and minimum number 
of connected pixels) were the same as used for the actual 
detections on the x 2 image. 

then perform aperture photometry on each separate fil- 
ter at the positions defined by the detection image. This 
procedure avoids the difficulties associated with merging 
separate single-band detection images (such as differences 
in object centroids between U— and I— band images, for 
example). This detection image is constructed using a \ 2 
image technique (Bzalay et al. 1999), expressed in equa- 
tion (g), where aj represents the background-subtracted 
pixel value in filter i, <ii the rms noise at this pixel and n 
is the number of filter, 



71 *■ ' 



(2) 



i=l 



This image has the advantage over other (more arbi- 
trary) combinations of images such as (V + I) in that 
it has a simple physical interpretation, namely that each 
pixel of this image represents the probability of detect- 
ing an object at that location. We compute aj and Oi 
for each pixel in each of the stacks using sextractor; 
this procedure allows us to correctly account for regions 
of varying signal-to-noise such as the overlap regions at 
the CCD boundaries. This resulting image is then used 
as input to sextractor as a detection image in the dual- 
image mode. We note that this method requires that both 
images are convolved to have the same full-width at half 
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Fig. 7. Comparison between CFDF and CFRS magni- 
tudes measured in 3" aperture diameter for the U-band 
(right panels) and I— band (left panels). The ~ 0.2 mag- 
nitude offset between CFRS and CFDF magnitudes for 
the 03hr field is a consequence of the significant galactic 
extinction in this field, which the CFRS magnitudes have 
not been corrected for. 



maximum and furthermore that they have a positional ac- 
curacy between filters of better than 1 pixel (as we have 



age. The filled shaded points show the median magnitude 
difference in half-magnitude intervals. Until Iab ~ 24, 
Iab (direct) - I AB (x 2 ) < 0.02; for 24 < I AB < 25.5, 
Iab (direct) — Iab(x 2 ) ^ 0.1. Beyond Iab ~ 24 magni- 
tudes computed using the direct image become systemat- 
ically brighter than the chisquared image, which is most 
likely a consequence of the more reliable object profile 
information contained in the chisquared image (which is 
comprised of effectively a sum of object fluxes over all fil- 
ters). In any case our galaxy colour measurements, which 
use aperture magnitudes, are unaffected by the applica- 
tion of this technique. 

Our final catalogues consist of matched V, I band cat- 
alogues for all fields. For fields 14hr, 22hr and 03hr we 
have additional B and U band imaging. Magnitudes in 
our catalogues are Kron ( Kron 1980| ) "total" magnitudes 
computed using the sextractor mag au to parameter. We 
have also carried out a comparison between these magni- 
tudes and those comp uted using the software employed in 
Le Fevre et al. (1986]) a nd the "Oxford" gala xy photome- 



demonfctrated m beet. 2~3 our internal positional accuracy 
is ~ 0.6 pixels, which meets tins objcctivej. We use an 
empirical approach to set the detection threshold in the 
chisquared image, similar to that employed in da Costa et 
al. (1998). Based on the numbers of objects detected in 
"blank" images (frames which have the same background 
noise as our real images) , the noise threshold is lowered in 
the chisquared image until the the number of additional 
sources detected is less than twice the number of sources 
detected in the blank images for the same change in the 
threshold. We emphasise however that the exact choice of 
the threshold is unimportant in this work as the range of 
variations considered in this procedure (~ 2. Oct) does not 
affect object detection even at the faintest magnitude limit 
where we carry out our scientific analysis (Iab ~ 25). 

3.2. Effect of x 2 -technique on galaxy magnitudes 

As object parameters crucial to galaxy photometry, such 
as half-light radius, are extracted from the detection image 
when using sextractor in dual-image mode (in addition 
to the normal (x,y) centroids) we wished to ensure that 
the use of the x 2 image did not bias our derived (total) 
magnitudes. In Fig. |§| we plot the difference in galaxy total 
magnitudes between the single band 03hr image and the 
dual-image mode method (% 2 image and stacked image) as 
a function of total AB magnitude in the single-band im- 



try package described in Metcalfe et al. (1991 ). We find no 
evidence of any systematic differences between these three 
softwares. Throughout this paper we measure colours in 
an aperture of 1.5" radius. 

We also perform star-galaxy separation using the rv, 
parameter from sextractor, which is carried out on the 
I— band catalogue. This parameter measures the radius 
which encloses half the object flux. Star-galaxy separation 
is not carried out faintwards of Iab ~ 21.5; in any case for 
high galactic latitude fields like ours galaxies outnu mber 
stars by a large fraction at these faint magnitudes (Rcid 



et al. 1996 ) . The bright limit of our catalogue (above which 
all galaxies and stars are saturated) is Iab ~ 18.5. 



4. Galaxy counts and colours 

4.1. Iab~ snd Bab — selected galaxy counts 

In Fig. |)| and Fig. [lC] we present I— and B— band galaxy 
counts as a function of AB magnitude derived from the 
CFDF sur vey (open thick circle s) compared with the com- 
pilation of Metcalfe et al. (2000| ) (small points). For the I— 
band we have four independent fields, or approximately 
~ 10 5 galaxies; for the B— band there are currently three 
fields with a total of ~ 0.75 x 10 5 sources. In both these 
filters we find an excellent agreement with the literature 
compilations extending over almost six orders of magni- 
tude. 

We find the fitted slope, a, for 20 < Iab < 24 to be 
0.3 5 ± 0.02, which agrees quite well with the value quoted 
in [Vletcalfe et al. (200C| ) of 0.33 for the slightly deeper 
limits of 21 < 1& C < 25. Note that all our /—data, from 
Iab ~ 18 to Iab ~ 24 (beyond which the effects of in- 
completeness begin to be important) is consistent with a 
constant slope in the number-magnitude relation. 



We also compare our counts to those of Postman et al 



1998| ) 



(filled squares). This dataset is a 16deg survey 
comprising 256 separate exposures of 900s each. This 
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Fig. 9. Tab— band galaxy counts for the four fields in 
the CFDF survey (open circles). The dashed line shows 
a least-squares fit to the bright data (22 < Iab < 24). 
The error-bar on each point corresponds to the field-to- 
ficld variance. The filled squares show galaxy counts from 
Postman et al. (1998|) ; the small points are from the com- 



pilation of Metcalfe et al. (2000). 



dataset is ~ 20% below the CFDF faintwards of Iab ~ 22 



and r emains 



until the effects of 



so mini Lnc eiiecis oi incompleteness ue- 
gin to ' dominate the CFDF counts. Fmthciinoic, in oiu 
source counts we see no evidence of the inflection seen at 
Iab ~ 22 in the Postman et al. 



data (Efstathiou et al. 1991; Roche et al. 1993; Braincrq 



et al. 1994; Hudon & Lilly 1996; Woods & Fahlman 1997 



McCracken et al. 200C). Here we will only briefly outline 



counts. 



the relevant equations. We have computed the two-point 



Fori the B-coimts we find a fitted slope iu the range 
20 < B A B < 24 of 0.47 ± 0.02, consistent with the value 



projected gal axy correlation function using the Landy & 



Szalay (1993 ) (LS) estimator, 



of 0.50 quoted in Brunncr et al. (1999) 



4.2. Galaxy colours 

In Fig. [ll] we show (V — I)ab colours measured as a func- 
tion of Iab total magnitude for galaxies in the CFDF llhr 
field (for clarity only 1/8 of all galaxies have been shown). 
The dashed line represents the colour incompleteness com- 
puted using the limiting magnitudes given in Table H. Our 
median (V — I)ab colours in half magnitude slices (filled 
circles) shows good agreement with colours derived from 
the NTT deep field (Arnouts et al. 1999b), shown as open 
squares, and this agreement continues to the completeness 
limit of our data. 



5. Measurements of u(6) in the CFDF fields 

5.1. Measuring u (6) 

There is an extensive literature on the measurement of 
the projected galaxy correlation function in deep imaging 



>{6) = 



DD - 2DR + RR 
RR ' 



(3) 



with the DD, DR and RR terms referring to the 
number of data-data, data-random and random-random 
galaxy pairs between 9 and 9 + 69. In this work we use 
logarithmically spaced bins with \og(S9) = 0.2. 

The fitted amplitudes quoted in this paper assume a 
power law slope for the galaxy correlation function, oj(6) = 
A^9~ ; however this amplitude must be adjusted for the 
"integral constraint" correction, arising from the need to 
estimate the mean density from the sample itself. This can 
be estimated as (Roche et al. 1993), 



C-± 



uj(9)dn 1 dQ 2 , 



(4) 



where ft is the area subtended by each of our survey 
fields. For the CFDF fields, we find C ~ AA^ by numeri- 
cal integration of Equation [| over our field geometry and 
assuming that galaxies closer than 1" cannot be distin- 



H. J. McCracken et al.: The Canada-France deep fields survey-I 



11 



1! 

-d 
\ 

60 

(0 






ao 3 

o 





1 1 1 1 I I 1 1 1 1 1 1 1 1 1 1 '1 1 ■ 








/ ■ ■ 


_ 


/ ■ ■ 




/■ .. . 




.7 • .! 




.. 4 . ■ 




° /- 




y 




. -■£ » ® $ 




. <? 




:¥ 




-9 


- 


.¥: 


- 


M- 




./' 




:& 


_ 


. >* O CFDF (this work) - 


- 


■ 2 


_ 


F 


- 


w 


1 


, :,/!' i i i 1 i i i 1 i i i 1 i i i 



18 20 22 24 

B.„ 



26 28 



Fig. 10. 

counts. 



As in Fig. g but for the Bab selected galaxy 



guished. In deriving the correlation amplitudes we then 
fit 



LO obs {0) = A^e- 



C. 



(5) 



where uj b s (8) is our observed measurement of u)(0). 
(We note however that the integral constraint correction 
only becomes important at larger scales (9 > 2') and that 
the large numbers of galaxies in our survey combined with 
the large field of view of each patch means we can fit for 
Aj(I') while neglecting the integral constraint, providing 
the range of the fit is restricted to scales where the inte- 
gral constraint is not i mportant). In our fitting procedure 
we use the method of Marquardt [1963 ) which takes into 
account that each measurement of ui(9) at each angular 
separation is not independent of the others. 

To minimise computational requirements (necessary 
given the large numbers of galaxies in our samples, typi- 
cally 30,000 galaxies per field) we use a sorted linked list 
method to compute the numbers of galaxy pairs in each 
angular bin (and furthermore we have verified that this 
method gives the same results as a direct-pair counting 
approach; it is, however ~ 2 — 3x faster). We use polyg- 
onal masks to blank out regions surrounding bright stars, 
large galaxies, satellite trails and cosmetic defects. We also 
discard data from CCD 8 because of its poor charge trans- 
fer properties. Our masking strategy is quite conservative, 
and (excluding CCD 8), around 15 — 20% of the total area 
is removed. Because the llhr and 22hr fields are composed 
of two separate stacks which are rotated relative to each 
other, are able to cover the entire 28' x 28' area for these 
fields. We have verified that the masks do not bias the 
determination of ui{6) function by applying them to ran- 
domly generated catalogues with the same number density 
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Fig. 11. Median (V — I)ab galaxy colour, measured as 
a function of Iab total magnitude, for galaxies in the 
CFDF- llhr field (for clarity only 1/8 of all galaxies are 
shown). The filled circles show the median galaxy colour 
in half magnitude intervals; the open squares show the me- 
dian (V — I)ab colou rs measured in the NTT deep field 
( Arnouts et al. 1999b). The dashed line shows the nominal 
colour incompleteness, based on the detection limits given 
in Table |. 



as the real catalogues and verifying that u>(6) is zero at all 
angular scales. We have also tested the effect of masks 
on clustered catalogues generated using the method of 
[Boncira fc Peebles (1978| ). These tests (which were carried 
out using the LS estimator) show that the application of 
masks does not change our fitted correlation amplitudes. 
Note that we do not apply a stellar dilution correction to 
our measured correlation amplitudes as some authors do, 
as the form of the stellar counts for high galactic latitude 
fields is not precisely known. However, this correction is 
unlikely to raise the amplitude of the correlation function 
by more than 10— 15%, given that stars outnumber galax- 
ies by at least an order of magnitude beyond Iab ~ 21.5. 

5.2. lo{9) for Iab magnitude-limited samples 

For each of our four fields, we divide our catalogue into 
magnitude limited samples. Each sample reaches one half- 
magnitude deeper than the previous, while the bright limit 
is kept at Iab = 18.5. (As the 14hr field is shallower 
than the other fields, we do not measure u>(9) on this field 
fainter than Iab — 24.0). These samples are extracted 
from the \ 2 catalogues prepared as outlined in Sect. |3.l| ; 
we have also performed the same computation on single- 
band catalogues (i.e, those computed without using an 
associated y^ image for detection) and find identical re- 
sults. Fig. [12| shows the logarithm of the amplitude of 
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Fig. 12. The logarithm of the average amplitude of the 
angular correlation function, ui{6), as a function of the log- 
arithm of the angular separation in degrees, for all fields 
in the CFDF. The error bar on each bin is computed from 
the field-to-field variance. Each of the five different sym- 
bols shows a different magnitude-limited sample; the solid 
lines shows the fitted correlation amplitude. This weighted 
fit is carried out neglecting the three innermost bins and 
assuming a power-law slope of —0.8 for uj(9) and a value 
of 4.2 for the integral constraint term C. Data from the 
14hr field is not included in the faintest slice. 



uj{&) averaged over the four fields of the CFDF as a func- 
tion of the logarithm of the angular separation in degrees, 
for a number of magnitude-limited samples. In the range 
-2.5 < log(0) < —1.3 (0.2' < 9 < 3') our measurements 
follow the expected power-law behaviour very well and at 
least until Iab ~ 24 there is no evidence of significant 
excess power on scales larger than the individual UH8K 
CCDs (~ 10'), in agreement with the evidence presented 



in Sect. 2.5 that the systematic photometric errors in our 



fields are negligible. 

At 9 < 0.2' we find at fainter magnitudes that u>{ff) de- 
viates from the expected power-law behaviour. We have 
attempted to investigate the origin of this effect (which is 
seen in all our fields) by carrying out an extensive set of 
simulations. In these simulations we generate a catalogue 
with ~ 40, 000 objects which have the same correlation 
amplitude as the Iab ~ 24 sample. Each of these objects 
is assigned a random magnitude in a specified interval 
and then added to the data frame. Object detection and 
photometry is carried out in a box extracted at this loca- 
tion. A new catalogue is constructed containing recovered 
total magnitudes for each object, and a magnitude cut 
is then applied to this sample. Objects which are lost in 
this process are typically those falling on or near bright 



Fig. 13. Simulated measurements of u)(6), (b,c), compared 
with CFDF measurements for slice 18.5 < Iab < 24.0, 
(a). In each case the simulated measurements have been 
normalised by the fitted amplitude for the 18.5 < Iab < 
24.0 slice, taking into account the integral constraint cor- 
rection and assuming a 5 — —0.8 power law. The sim- 
ulations (b,c) contain objects in the magnitude range 
23.5 < I A b < 24.5 and 24.5 < I AB < 25.5 respectively. 
Each plot has been offset by an arbitrary amount x where 
x = 0,1.5,3.0. 

stars or galaxies, or those whose recovered total magnitude 
falls outside our magnitude cut. Masks are applied to this 
catalogue and u>(9) computed using the same procedures 
used for the real data. This procedure is then repeated 
for progressively fainter magnitude slices. (In construct- 
ing the simulated catalogue, two simplifying assumptions 
were made, firstly that the input magnitude distribution of 
objects is flat, and secondly we use objects with Gaussian 
point-spread functions.) 

In Fig. [13| we show the results of two of these simula- 
tions, labelled (b) and (c). The results displayed for (a) 
show the measured correlations for the 18.5 < Iab < 24.0 
magnitude slice. In order to display deviations from the 
power-law behaviour, each slice has been normalised by 
the fitted amplitude for the 18.5 < Iab < 24.0 slice, 
taking into account the integral constraint correction and 
assuming a 5 = —0.8 power law. These simulations can- 
not reproduce the depression in the correlation function 
found at scales log(0) ~ —2.5. We have also investigated 
if the depression is caused by an excess of objects around 
bright stars, and have found no evidence of such an ex- 
cess. An important aspect of these simulations, however, 
is that they confirm that excess power on large scales only 
becomes important for the CFDF data for the faintest 
magnitude ranges, 24.5 < Iab < 25.5. 
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Fig. 14. The logarithm of the fitted amplitude of the angular correlation uj{9) at 1' as a function of sample median 
Iab magnitude (filled circles). Fits were computed assuming a power law of slope 8 — —0.8 for lo(9) and an integral 
constraint term calculated as described in the text. Error bars on the CFDF points are computed from the analytic 



expression of Bernstein (1994). The model curves are taken from McCracken et al. (2000) and show three different 
clustering evolution scenarios which are detailed in the text. 



In Fig. |lj we plot the fitted amplitude of uj(9) at 
one arcminute, A w (l'), averaged over all our fields, (filled 
circles) as a function of the sample median AB magni- 
tude {I AB-med)- We have estimated the error bars on 
A u , 8A U , using an analytic approximation introduced by 



Bernstein (1994 ) and further developed by Bzapudi fc 
Colom fi (1996 ); see also Arnouts et al., in preparation 



for a detailed exp lanation o f the application of this ap- 
proximation. The Bernstein estimator relies on a knowl- 
edge of the higher-order moments of the galaxy correla- 



tion function, S3 and S4; we have estimated these quanti- 
ties directly from the CFDF dataset and they will be pre- 
sented in a future work (Colombi et al., in preparation). 
At bright magnitudes (I AB-med ~ 18), SA^, is dominated 
by the Poisson error component; however, faintwards of 
lAB-med ~ 22, the main component of SA^ in our field 
consists of cosmic variance (or "finite volume") effects. 
Our analytic estimates of this effect indicate that the er- 
rors estimated empirically from the field-to-field variance 
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of our four fields may underestimate the total error at 
these magnitudes by around ~ 20%. 

Our fitted amplitudes were computed assuming a —0.8 
slope for w{9) (to allow a comparison with the literature), 
and over the range 0.2' < 6 < 3'. We have also attempted 
to fit for A^ only at larger separations (~ 3'), but our 
survey fields are not large enough to detect any scale- 
dependence in the A^ — iAB-med relation. Additionally, 
Fig. R4 shows a compilation of recent measurements of 
Au from the literature. These literature measurements all 
assume a fixed slope of S = —0.8, with the exception of the 



Postman et al. survey, in which the slope was allowed to 



vary with the fit. To allow a fair comparison with the other 



authors, and with our work, the Postman et al. points are 
not plotted faintwards of IAB-med ~ 21, where their fitted 
slopes begin to differ significantly from —0.8. 

At bright magnitudes (19 < IAB-med < 22) we find 
that our observations are compatible with almost all the 
data from the literature compilation. At fainter magni- 
tude ranges, (22 < IAB-med < 24), our observations 
clearly favour a low amplitude for A u , At Iab ~ 24 they 
are least a factor of ten below the value of A^ found by 
Braincrd fc Smail (199S ). This work covers a small area 



50 arcmin ) and consists of indivi dual unconnected 



pointings. As has been suggested before ( [McCracken et al 
2000|)t]he discrepancy may be due to field-to-field variance 
in the galaxy clustering signal. To provide a more direct 
answer this question, we extract 200 regions of 50 arcmin 2 
from the 22hr and llhr fields (these two have complete 
coverage as a consequence of bonette rotations). On each 
of these sub-regions we measure A^ as for the full sample. 
In Fig. [15] we show the histogram of the fitted values for 
the llhr and 22hr fields; from this we estimate that a ±3er 
error bar indicates a dispersion of xlO in A u . It is clear 
that the error bars on Brainerd & Smail measurement un- 



derestimate the true error by a large amount. 

We have also determined A w (l') in one-magnitude 
slices, for instance, 19.5 < I A b < 20.5, 20.5 < Iab < 21.5 
to the limit of our survey, as an additional check. These 
measurements are considerably more noisy than the in- 
tegrated measurements presented above because of the 
smaller numbers of galaxies in each slice. However, we find 
that the derived IAB-med ~ log Aa (lQ relation is very sim- 
ilar to what is presented in Fig. [L4[ 

5.3. Dependence of slope on magnitude 

The large numbers of galaxies in our survey allows us to 
make a direct measurement of S, the slope of u>(6) as a 
function of sample limiting magnitude. Several attempts 
at this measurement have been carried out but with gener- 



ally inconclusive results. Neuschaefer & Windhorst (1995) 



find 5 = — 0.5 at g < 25 based on two independent fields 
each the size of one of our CFDF fields. Shallower wide- 



angle surveys like th ose of Roche fc Eales (1999 ) and 
Cabanac et al. (200C ) find no evidence for deviation from 
a slope of ~ 0.8. One difficulty in this measurement is that 




Fig. 15. The logarithm of the fitted correlation amplitude 
at 1', log(Aui) measured on 200 7' x 7' sub-areas extracted 
from the llhr and 22hr fields. Galaxies extracted have 
magnitudes in the range 18.5 < Iab < 25.5. The dot- 
ted lines represent ±3<r confidence limits about the mean 
value, shown by the solid line. The arrow shows approxi- 
mately the A u (l') measured by Braincrd & Smail (1998). 



excess power on large scales, such as can be produced by 
zero-point variations, can produce an artificially shallow 
slope. However, we have already demonstrated that our 
magnitude zero-point errors are small across our mosaics 
(Fig. |]) and that differential incompleteness only becomes 
significant for measurements of u>(9) in the CFDF fields at 
very faint magnitudes (Fig. |l3|). Nevertheless, we adopt a 
cautious approach and use two different methods to mea- 
sure 5. 

Firstly we compute \ 2 contours from the average corre- 
lation per angular bin over all the fields (Fig. |lj), shown 
in Fig. nq. In our second method we compute the slope 
independently for each field and measure the field-to-field 
standard deviation of this fit, which is shown in Fig. n7L 
together with a comparison with the fitted slopes from 



Postman et al. (1998), computed for a fitting range of 



0.5' < 9 < 5'. An additional complication is that in both 
cases our slope-fitting procedure requires an estimation 
of the integral constraint (Equation ra), which in turn de- 
pends on the slope. To limit these difficulties, we perform 
the fit only in the range —0.2' < 9 < 1.9' where the effects 
of the integral constraint are expected to be negligible. 

Both methods indicate that at bright magnitudes, 
18.5 < Iab < 22.0, 6 ~ —0.8; at fainter magnitudes we 
detect a slight flattening of the slope, with 6 ~ —0.6. The 
error bars in Fig. |l7| show the error in 5 for a given value of 
A u ; from Fig. [U] we see, however, that a slope of S = —0.8 
is still within 3<r of our best fit for all magnitude ranges. 
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Fig. 16. x 2 contours showing best-fitting amplitudes and 
slopes (plus symbols) for the four faintest CFDF samples, 
from right to left, 18.5 < Iab < 22, 23, 24, 25. The three 
contours show the lcr (thick contour), 2<r and 3cr confi- 
dence levels. 
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Fig. 17. The fitted slope of ui(9) as a function of sample 
median magnitude (open circles) for galaxies with 18.5 < 
Iab < 22, 23, 24, 25 averaged over the CFDF fields. Points 
are also shown from Postman et al 1998 (filled circles). 
Error bars on the CFDF points are computed from the 
variance over our four survey fields. 

5.4. Galaxy clustering as a function of colour 

For all galaxies in our sample we can measure (V — I)ab 
colours and we can use this to select subsamples by colour. 
In Fig. [L8| we show log A„(l') as a function of (V — I)ab 
colour for galaxies with 18.5 < Iab < 24.0 (open pen- 
tagons) and 18.5 < Iab < 23.0 (asterisks). For each of 
the four fields, we divide our sample into twelve equally 
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Fig. 18. The logarithm of the fitted amplitude of the an- 
gular correlation u>(6) at 1' as a function of sample me- 
dian Iab magnitude (as before, fits assume 5 = —0.8). 
\og(A UJ )(l') was computed for 12 equally spaced bins of 
0.25 in (V — I)ab- Error bars for each colour slice are 
computed from the field-to-field variance. Two magnitude 
ranges are shown, 18.5 < Iab < 24.0 (open pentagons) 
and 18.5 < Iab < 23.0 (asterisks), together with the fit- 
ted amplitude for the full field samples, dotted and dashed 
lines respectively. 



spaced bins in colour each of width 0.25 in (V-I)ab- The 
error bars in logA^l') were computed from the field-to- 
ficld variance. Also shown as the dotted and dashed lines 
is the amplitude of log A U1 {Y) for the full-field sample for 
these two slices. 

We first note that both these cuts are relatively bright; 
at Iab ~ 24 our catalogues are expected to be ~ 100% 
complete (Figs. [| ||). Additionally, as shown in Fig. O, 
the effect of colour incompleteness should be minimal, al- 
though at Iab ~ 24 we may begin to lose some extremely 
red objects from our sample. 

Our measurements clearly show that redder objects 
are more strongly clustered, as has been widely reported 
for local populations ( Loveday et al. 1995| ). We find that 
objects with (V — I)ab ~ 3 have a clustering ampli- 
tude at least a factor of ten higher than the full field 
population, shown as the dotted and dashed lines in 
Fig. 111. Interestingly, we also find that the bluest objects 
in our survey, those with (V — I)ab ~ 0, have a clus- 
tering amplitudes marginally higher than the full sample. 
Furthermore, we find that none of our colour selected sam- 
ples have clustering amplitudes below the full-field value. 
We discuss the implications of these results in the follow- 
ing section. 
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6. Discussion and comparison with model 
predictions 

6.1. Modelling A w 

For r < 20h~ 1 Mpc, the spatial correlation function 
£(r) can be approximated by £(r) = (ro/r) 7 , where, 
from the results of local redshift surveys, 7 ~ —1.8 and 
r ^ASh-i-Mpc (|Groth fc Peebles 1977J |Davis fc Peebles 
1985; Maddox et al. 1990J ). One way to produce a pre- 



diction for the variation of A w with sample limiting mag- 
nitude is to assume a functional form for the growth of 
clustering £(r, z), normally written as 



£(r, *) = /*(*) (7) 



where 



h[z) = (l + z)" (3+e) . 



(6) 



(7) 



This functional form is then integrated over redshift 
space using the relativistic version of Limber's equation, 



(Phillipps et al. 1978; Groth fc Peebles 1977; Limber 



1953), assuming that for $ < 1 (Efstathiou et al. 1991), 



u(9) 



-r[( 7 -l)/2] A 



9~ 



_1 '0' 



(8) 



r(7/2) 

where T is the incomplete gamma function, 9 is the angu- 
lar separation and A is given by 

A = 



with 



9{z) 



00 

9(z) 


/dN\ 2 
\~dz~) 


dzj 


h{z) 


(dr(z)Y 1 


d\-\; 


I dz 





dN 

dz 



i 2 



dz 



(9) 



(10) 



where d^(z) is the angular diameter distance and dr{z)/ dz 
is the derivative of the proper distance. In this simple 
scenario, three cases are of interest: clustering fixed in 
proper coordinates, in which e = 0.0; clustering fixed in 
co- moving coordinates which gives e = —1.2. Finally, the 
predictions of linear theory give e = 0.8. 

Many papers have investigated the scaling of A u with 
magnitude using the approach detailed above (see, for ex- 
ample, Efstathiou et al. (1991)). To interpret measure- 



ments of A U) using these models, however, involves making 
at least two critical assumptions: firstly, the form of the 
redshift distribution dN / dz for the faint galaxy popula- 
tion and its evolution as a function of limiting magnitude; 
and secondly how £(r, z) scales with redshift (Equation Q). 
In the following section we examine these two assumptions 
in turn. (Predicted correlation amplitudes using this for- 
malism are also sensitive to the underlying cosmology, as 
the size of the volume element at a given redshift is much 
lower for an Einstein-deSitter cosmology than for a low-17 
universe. However, to the median redshift of our survey 
the difference in model predictions between open and flat- 
Lambda cosmologies small. In this paper we assume that 
fl]\i ~ 0.1, in agreement with recent observational evi- 
dence.) 



6.2. Validity of model assumptions 

Our dN J dz is derived from luminosity evolution models 
which are described fully in Metcalfe et al. (2000| ). Starting 
with the observed local galaxy luminosity function and 
assuming a star-formation history for each galaxy type 
these models are able to reproduce the observed num- 
bers counts, colours and redshift distributions of the faint 
galaxy population to the limits of the current observations 
(Metcalfe ct al. 2000). However, at Iab ~ 24 we may now 



directly test these model redshift predictions against spec- 
troscopic measurements made in the Hubble Deep Field 
(HDF) ( pohen ct al."2000[ |Williams ct al. 1996| ). In Fig. [TJ 
(upper panel) we show the spectroscopic redshift distri- 
bution for 120 galaxies in the HDF-North (with 16 non- 
detections represented as the open box) compared to the 
predictions of our fl^j =0.1, ft\ = 0.0 PLE model distri- 
bution (solid line). In the lower panel we show the relation- 
ship between median redshift z me d and Iab limiting mag- 
nitude as predicted by this model (note that we compute 
our median redshift in each case by considering all galax- 
ies brighter than the abscissa magnitude). Additionally, 
we show measurements from several other Iab~ magni- 
tude limited redshift surveys, including the HDF sample 
used in the top panel. We also show the median redshift 
derived from photometric redshifts for two samples lim- 
ited at I A b < 25 for the HDF N/S, kindly supplied to 
us by S. Arnouts. The Poissonian error bars computed 
from all surveys are smaller than the symbols and are not 
plotted. The true field-to-field variance may of course be 
much larger, but the fact that we measure the same me- 
dian redshift for both HDF-N/S catalogues suggest that 
it is not. 

Despite these qualifying remarks, we conclude that 
our luminosity evolution models provide an acceptable 
fit to the observed redshift distributions at least to the 
depth to which we measure galaxy clustering in the CFDF 
photometric catalogues (lAB-med ~ 24, or equivalently 
Iab ~ 24.5). They are able to reproduce both the trend 
of z me d with Iab and the dispersion in redshift at a 
given magnitude slice. Furthermore, as demonstrated in 
Metcalfe et al., they also correctly predict the numbers of 
2 < z < 3 galaxies. We therefore conclude that our mod- 
elling of dN / dz is not a major source of uncertainty in our 
prediction of A^ . 

Our second assumption, that the growth of galaxy clus- 
tering can be expressed as in Equation (m), is more prob- 
lematic. Clustering measurements of Lyman-break galax- 
ies (Adelberger et al. 1998; Giavalisco et al. 1998| ), have 
already indicated that the "epsilon" formalism does not 
provide an acceptable fit to the observations . Similar re- 
sults have also been found for measurements of ro(z) in 
the HDF-North ( |Arnouts ct al. 1999a| ). Can our cluster- 
ing measurements in the CFDF be successfully matched 
by this model? In Fig. |20| we show our measurements 
of A U (V) compared to prediction of our models for e = 
— 1.2,0.0,1.0,2.0 (long dashed, solid, dotted and dashed 
lines respectively), assuming tq = 4.3/i _1 Mpc, Q,m = 
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Fig. 19. Upper panel: the redshift distribution of faint 
galaxies in the Hubble Deep Field North from the compi- 
lation of Cohen et al. (2000) (histogram), compared to the 
predictions of our model (smooth solid line). The open box 
represents the number of galaxies observed for which no 
redshift could be determined. Lower panel: the relation- 
ship between median redshift z me< i predicted by our model 
(solid line) compared to the CFRS (five-pointed stars), 
spectroscopic redshifts in the HDF-N (cross) and photo- 
metric redshifts from the HDF N/S (asterisk). Poissonian 
error bars are smaller than the symbols in all cases and 
are not plotted. 



0.1, fl A = 0.0 and S = -0.8. (We note that cluster- 
ing predictions for an Q — A cosmology are very similar 
to zero- A cosmology). For clarity we omit the literature 
compilation shown previously. In the magnitude range 
18.5 < lAB-med < 22, we see that our observations are 
consistent with e ~ 0. However, faintwards of Iab ~ 22, 
our observed clustering amplitudes decline more rapidly 
than the model predictions. By 23.0 < Iab < 24.0 our 
observations are consistent with e > 1. From Fig. Effl it 
is clear that the e ~ model cannot match simulta- 
neously both bright and faint observations in the range 
18.5 < lAB-med < 24.0. Furthermore, rapid growth of 
clustering for the entire sample (e ~ 2) is marginally ex- 
cluded because it produces correlations which are already 
too low by lAB-med ~ 22 to match our observations. 
Furthermore, allowing ro (i.e., ro(z = 0)) to vary merely 
changes the normalisation at Iab ~ 18.5 (which is already 
in agreement with our observations) b ut n ot the slope of 
the A^(l') — lAB-med relation. In Sect. ^J we investigate 
the reason for this discrepancy in more detail. 



6.3. Colour selected galaxy clustering 

In Fig. [l8| we clearly see the dependence of log(A w ) on 
(V — I)ab colour. To interpret this result, we first note 
that the dependence of (V — I)ab colour on redshift and 
morphological type is well established, thanks to exten- 



sive spectroscopic surveys (Lin et al. 1999; Cowie et al 



1996| ; |Crampton et al. 1995J ; prilly et al. 1995J ). In par- 



ticular, Wilson et al. (2001| ), using a large spectroscopic 
sample, demonstrated that objects with (V — I)ab ~ 3 
are predominately massive elliptical galaxies at z ~ 0.8. 
Furthermore, clustering amplitudes have recently been 
measured for objects selected to have extremely red 
colours in optical-infrared bandpasses ( Daddi et al. 2000 ). 
These objects have clustering amplitudes ~ 10 x higher 
than the full field population. It is probable that these 
objects are closely related to our (V — I)ab ~ 3 sam- 
ple; for galaxies with (V — I)ab > 3 at Iab < 23.0 , 
we find ~ 0.3 N ga/ arcmin -2 . In [Daddi et al. (2000| ) 
using a selection of (R — K s ) > 5.0 and the slightly 
brighter limit of K s < 19.0, they find a surface density 
of 0.5 N„ a ; arcmin~ 2 . The difference between our full field 



23 and the clus- 
3 is approx- 



clustering amplitude at 18.0 < Iab < 
tering of objects selected with (V — I)ab 
imately the same as the difference found by Daddi et al 
between their i-sT s -selected sample and those of their ex- 
tremely red objects. 

Intriguingly, at the blue end of our selection, (V — 
I)ab ~ 0, we also find a higher clustering amplitude than 
the full field sample, although the error bars are large due 
to field-to-field variations (at the ~ 0.1 magnitude level) in 
galaxy colours and the small numbers of objects involved. 
We have repeated our measurement of A u (\ r ) using an 
integrated selection (i.e., considering only objects redder 
or bluer than a specified colour cut) and find a similar 
effect. There is some evidence for this effect in the liter- 
ature: working with photographic data, and considering 
objects in a somewhat brighter blue selected magnitude 
cut, 20 < Bj < 23.5, Landy et al. (1996) also found an 



enhancement of at least ~ 10 for the clustering amplitudes 
of the bluest objects in (U — Rf). 

Although Lyman-break galaxies are expected to be flat 
spectrum objects and therefore have (V — I)ab colours of 
~ their surface densities are not large enough to pro- 
duce the effect seen in Fig. hS. The most likely explana- 
tion of this result is that these objects constitute a low- 
redshift population whose higher correlation amplitudes 
are a consequence of the lack of projection effects which 
dilute the measured A^s. Some evidence for this can been 
seen in Fig. 5. of Crampton et al. (1995J ); all objects with 
(V — I)ab ~ are at z < 0.3. We also note that our red 
and blue samples have very low cross-correlation ampli- 
tudes, which supports the notion that the objects in our 
survey with (V— I)ab ~ 3 and (V — I)ab ~ are separate 
populations at different redshifts. 
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mass, show clearly that more luminous galaxies have clus- 
tering amplitudes very different from less luminous ones 



(Benson et al. 2001; Baugh et al. 1999; Kauffmann et al 



Furthermore, it is now reasonably well established 
from observations that locally the galaxy correlation 
length tq depends on morphological type and colour 



(Tucker et al. 1997; Loveday et al. 1995; Davis & Gcllci 



1976), and some evidence exists for a direct dep endence 



between lu minosity and clustering amplitude (Bcnoist 
1996). There are also indications that these trends 



Fig. 20. The evolution of log(Aj(l')) (assuming S — 
-0.8, fi^ = 0.1 and r = 4.3/i,- 1 Mpc) for e = 
— 1.2,0.0,1.0,2.0 (long dashed, solid, dotted, and dashed 
lines respectively). The filled points show the CFDF full 
sample. This plot is similar to Fig., [lj except the literature 
compilation has been omitted for clarity. 



6.4. Biasing and the growth of galaxy clustering 

Given that the redshift distribution used in our models is 
in agreement with our observations, then it is clear that 
the discrepancy evident in Fig. EG between model predic- 
tions and observations must be a result of evolution in the 
intrinsic properties of the galaxy population. Our simple 
model does not take this into account. As a first step to- 
wards a more realistic description of the data, we may 
try changing the form of the Tq 



continue to higher redshifts ( Carlberget al. 2000 ; Le Fevre 
et al. 1996 ). Furthermore, in our dataset, the median field 
galaxy (V — I)ab colour changes by ~ 0.4 mag in the 
range 22 < I AB < 24 (Fig. |l|); from Fig. |l| we see that 
changes in colour of this magnitude cannot produce the 
changes in amplitude of .4^(1') seen in the data. For this 
reason we suggest that the rapid decline in A^(l') in the 
range 22 < lAB-med < 24 is a consequence of luminosity- 
dependent clustering segregation. 

Extensive imaging and spectroscopic observations have 
demonstrated that, for any magnitude limited sample, as 
we probe to fainter magnitudes, the mean intrinsic lumi- 
nosity of the field galaxy population becomes progressively 
fainter. This is illustrated in Fig. E^ where we show the 
absolute luminosity as a function of apparent magnitude 
for galaxies in the CFRS survey (Lilly et al. 1995) (open 
circles) and for galaxies in the HDF-North computed us- 
ing photometric redshifts from the photometric catalogue 



z relation: McCracken 
et al. ( p000| ), considered such a modification by adopting 
the form of £(r, z) derived from dark matter haloes with 
Vmax > 120 kms -1 i dentified in a large, high-r esolution 
N— body simulation ( Kravtsov fc Klypin 1999 ). This is 



shown as the solid line in Fig. |l4|. However, because the 
form of this relationship is very similar to the traditional 
"epsilon- model" in the range < z < 1, and because there 
are few z > 3 galaxies in samples limited at Iab < 25, the 



of Fcrnandez-Soto et al. (199E). For both these catalogues 
we also show the median absolute magnitudes computed 
in half- magnitude intervals of apparent magnitudes (open 
and filled circles for the CFRS and HDF respectively). 
We also include an estimate of the median differential lu- 
minosities dL/dm from our luminosity evolution models 
(solid line). We see that between Iab ~ 20 and Iab ~ 22, 
within the CFRS sample, the median galaxy luminosity 
declines by ~ 0.5 magnitudes. In the range 22 < Iab < 24 
we measure a decline of a further magnitude, although 
the uncertainties in the absolute luminosities computed 
from photometric redshifts in the HDF is probably at least 
^0.5 magnitudes. The decline in model luminosities seen 
in the range 18.5 < Iab < 22.5 are a consequence of the 
steep faint-end slope we adopt for the spiral galaxy lumi- 
nosity function. These faint galaxies at I ~ 22 are pre- 
dominately late-type galaxies, as has been demonstrated 
by spectral and morphological classification (Brinchmann 



differences in predicted amplitudes between this and the et al. 1995; Driver et al. 1998) 



conventional formalism are small in the magnitude ranges 
we consider. 

The basic reason why the "e" -models fail to repro- 
duce the clustering of Lyman-break galaxies and the ob- 
served form of r (z) at high redshift is that they implic- 
itly ignore the existence of bias and that how galaxies 
trace the underlying d ark matter depends on the mass o f 
the dark matter halo ( Kaiser 1984 ; Bardeen et al. 1986| ). 
Measurements of galaxy clustering in semi-analytic mod- 
els, which provide a prescription for how galaxies trace 



7. Conclusions 

In this paper we have introduced a new wide-field deep 
imaging survey, the Canada-France deep fields project. 
Our survey covers a total area of ~ 1 deg 2 in four separate 
contiguous 0.25 deg 2 fields. We have demonstrated that 
our reduction procedures are robust. Our internal astro- 
metric errors are ~ 0.06" and our systematic photometric 
errors across each 0.25 deg 2 field are < 0.04 magnitudes. 
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Fig. 21. Absolute J^s magnitudes (M/) and median ab- 
solute luminosities as a function of Jab magnitudes for 
galaxies selected from the CFRS (open circles, stars) and 
from the HDF-N (small and large filled circles). Model 
predictions of our luminosity evolution model are shown 
as the solid line (In all cases luminosities are computed 
assuming h — 0.5 and an open Qm = 0.1 cosmology). 



In each of our four fields our galaxy catalogue is 
- 100% complete for 18.5 < I AB < 24.0. We find that 
in this range the Iab number-magnitude relation is well 
fitted by a line of constant slope dlog(N)/dm ~ 0.35±0.02 
in agreement with the literature compilations. Our com- 
pleteness drops to 50% at Iab ~ 25.5. 

Using this survey we have measured the projected-two 
point correlation function u>[&) for a sample of 100,000 
galaxies as a function of sample median magnitude, 
lAB-med, angular separation 6 and (V — I)ab colour to 
iAB-med ~ 25.0. At lAB-med ~ 24 we measure the am- 
plitude of w{6) at 1', Aj(I') to an accuracy of 30%. Our 
conclusions are as follows: 

1. We find that in the range 19.5 < lAB-med < 24.0, 
A^, declines monotonically with sample limiting magni- 
tude and that throughout this range, u)(9) is well matched 
with a power-law of slope 6 for 0.2' < 8 < 3.0'. At 
bright magnitudes, S ~ —0.8; at lAB-med ~ 23, we find 
8 ~ 0.6, although our observations are still compatible 
with 5 ~ —0.8 at a 3a confidence level. 

2. We find a clear dependence of A u (l') on (V — I)ab 
colour for galaxies selected in two magnitude ranges, 
18.5 < Iab < 24.0 and 18.5 < I AB < 23.0. Galaxies with 
(V — I)ab ~ 3 have clustering amplitudes ~ 10 times 
higher than the full field population. These objects are 
most probably evolved ellipticals at z ~ 1. We also find 
some evidence (at the ~ la — 2a level) for slightly higher 
clustering amplitudes for the blue (V — I)ab ~ objects 
in our sample. 



3. We discuss model predictions and current spectro- 
scopic determinations of the redshift distribution dN/dz 
for the faint field galaxy population. We conclude that 
for 19.5 < lAB-med < 24.0, dN/dz is now well deter- 
mined. Using these predictions we find that for low f2 cos- 
mologies, and assuming a local galaxy correlation length 
ro = 4.3ft.~ 1 Mpc, the growth of galaxy clustering, param- 
eterised by e, is consistent with e ~ for galaxies in the 
magnitude range 19.5 < lAB-med < 22.0. 

4. However, in the magnitude range 22.0 < lAB-med < 
24.0, our observations are consistent with e > 1. Models 
with e ~ cannot match simultaneously measurements of 
Aj(I') at bright (lAB-med ~ 19) and faint {lAB-med ~ 
24) magnitudes. 

5. We demonstrate that one simple interpretation of 
this result is that by lAB-med ~ 24 our sample is dom- 
inated by intrinsically faint (Mi ~ —20) galaxies which 
have considerably weaker correlation lengths (ro ~ 2 — 
3/i -1 Mpc) than the local galaxy population. 
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